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Abstract 

For any prime power q and any dimension s > 1, we present a 
construction of (i, s)-sequences in base q with finite-row generating 
matrices such that, for fixed q, the quality parameter t is asymptoti- 
cally optimal as a function of s as s — > oo. This is the first construc- 
tion of (t, s)-sequences that yields finite-row generating matrices and 
asymptotically optimal quality parameters at the same time. The con- 
struction is based on global function fields. We put the construction 
into the framework of (u, e, s)-sequences that was recently introduced 
by Tezuka. In this way we obtain in many cases better discrepancy 
bounds for the constructed sequences than by previous methods for 
bounding the discrepancy. 

1 Introduction and basic definitions 

Constructing sequences with good equidistribution properties is an impor- 
tant problem in number theory and has applications to quasi-Monte Carlo 
methods in numerical analysis (see, e.g., [31 EE]). In this context, the star 
discrepancy appears as an important measure of uniform distribution. For a 
given dimension s > 1, let J be a subinterval of [0, l] s and let x Q , . . . , x^-i be 
N points in [0, l] s (we speak also of a point set V of N points in [0, 1} S ). We 
write A(J; V) = A(J; x , . . . , a; A r_ 1 ) for the number of integers < n < N — 1 
for which x n G J. Then the star discrepancy of the point set V consisting of 
the points Xq, . . . , x^^i is denned by 



D* N (P) = D* N (x , . . . , ccjv-i) = sup 

j 
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where the supremum is extended over all subintervals J of [0, l] s with one ver- 
tex at the origin. For a sequence S of points x , x%, . . . in [0, l] s , the star dis- 
crepancy of the first N terms of S is defined as D^(S) = D^(x , . . . , x^-i)- 
We say that S is a low- discrepancy sequence if 

D* N (S) = 0(N- 1 (log N) s ) foralliV>2, 

where the implied constant does not depend on N. This is conjectured to be 
the least possible order of magnitude in N that can be obtained for the star 
discrepancy of a sequence of points in [0, 1} S . 

One of the most powerful methods for constructing low-discrepancy se- 
quences is built on the theory of (t, s)-sequences and (t, m, s)-nets which was 
developed by Niederreiter [13] on the basis of earlier work by SoboF [22J and 
Faure [3J. The reader is referred to, e.g., the monographs [3 J and [16] for the 
general background on this theory. In the following, we give a short descrip- 
tion of the basic notions. For integers b > 2 and < t < m, a (t, m, s)-net in 
base b is a point set of b m points in [0, l) s such that every elementary interval 
J Q [0, l) s in base b with volume b l ~ m contains exactly b l points of the point 
set. By an elementary interval in base b we mean an interval of the form 

s 

W[ ai b- d %(a t + l)b^) 
i=i 

with integers di > and < a, < b di for 1 < i < s. A sequence Xq, Xi, . . . 
of points in [0, l] s is called a (t, s) -sequence in base b, where b > 2 and t > 
are integers, if for all integers k > and m > t the points [jE n ]b,m with 
kb m < n < (k + l)b m form a (t, m, s)-net in base b. Here [aJ n ]b,m denotes the 
coordinatewise m-digit truncation in base b of the point x n . (For a detailed 
description we refer the reader to, e.g., [24], [251 Section 6.1.1], and [271 
Section 2] ; the important role of the truncation was emphasized again in the 
recent paper [5].) In this concept the integer t serves as a quality parameter, 
and generally speaking the smaller t the more uniform the distribution of the 
sequence. 

Most of the known (t, s)-sequences are obtained by the so-called digital 
method which was developed by Niederreiter [13] and constructs a sequence 
as follows. Choose a dimension s G N, a finite field ¥ q with cardinality q, 
and put Z q — {0, 1, . . . , q — 1} C Z. Note that q must be a prime power here. 
Choose 

(i) bijections ip r : Z q — > ¥ q for all integers r > 0, satisfying ip r (0) = for 
all sufficiently large r; 

(ii) elements c^ r e ¥ q for 1 < i < s, j > 1, and r > 0; 
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(iii) bijections Ajj : ¥ q — > Z q for 1 < i < s and j > 1. 

For the construction of the sequence we make use of the notion of the gen- 
erating matrices := {c^\)j>i, r >o £ F^ xN ° for 1 < z < s. If these matrices 
. . . , satisfy the property that each row of each matrix contains only 
finitely many nonzero entries, then we speak of finite-row generating matrices. 
The ith coordinate ii'' of the nth point x n = (xn \ • • • , %') of the sequence 
is computed as follows. Given an integer n > 0, let n = Yl^Lo z r( n )q r be 
the digit expansion of n in base q, with all z r (n) G Z q and z r (n) = for all 
sufficiently large r. Then for 1 < i < s we form the matrix- vector product 
over F q given by 







/v (i) \ 


C (i) . 




(i) 

— '■ Vn,2 




V / 


\ J 



Finally, we put 

oo 
3=1 

The distribution of the sequence Xq, x±, . . . is mainly determined by the 
rank structure of the generating matrices. It is well known that the digital 
method generates a (t, s)-sequence in base q if for every integer m > t and all 
nonnegative integers d±, . . . ,d 8 such that di + - ■ ■ + d s = m— t, the (m-t)xm 
matrix over ¥ q formed by the row vectors 

with 1 < j < di and 1 < i < s has rank m — t (see [21 Section 4.4] and [TBI 
Section 4.3]). 

Standard (t, s)-sequences obtained by the digital method include the 
Sobol' sequences [22], Faure sequences [1], Niederreiter sequences [Hj, gener- 
alized Niederreiter sequences [24] , and Niederreiter-Xing sequences [T8], [27] . 
Summaries of the constructions of these sequences can be found in [3], Chap- 
ter 8] and [UJ. 

The digital method can be applied also for the construction of (t, m, s)- 
nets in base q. In this case, the generating matrices . . . , C^ s ' are m x m 
matrices over ¥ q (see [3l Section 4.4.1]). 

Tezuka [26J recently pointed out a deeper regularity in the generating 
matrices of the generalized Niederreiter sequences and introduced so-called 
(w, m, e, s)-nets and (w, e, s)-sequences. He also established a general discrep- 
ancy bound for (u, e, s)-sequences and in this way he improved the constants 
in the discrepancy bounds known for generalized Niederreiter sequences. 
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In this paper we fulfill several objectives. In Section [5] we give a revised 
version of Tezuka's definition of (u, m, e, s)-nets which allows us to determine 
a corresponding quality parameter t for a (u, e, s)-sequence in the language 
of (t, s)-sequences. In Section [3] we introduce a new construction principle for 
generating matrices of digital (w, e, s)-sequences using global function fields. 
From the viewpoint of (t, s)-sequences, these new sequences have the same 
quality parameter t as the Niederreiter-Xing sequences, but the construction 
is simpler. Moreover, a refined construction in Section H] yields (u, e, s)- 
sequences and (t, s)-sequences with the same parameters as in Section [3], but 
with finite-row generating matrices. Finally, in Section [5] we discuss upper 
bounds on the star discrepancy for the sequences constructed in this paper. 

2 A revised definition of (u, e, s)-sequences 

We follow up on the recent work of Tezuka [26] on nets and (t, s)-sequences, 
but for several reasons (for instance, in order to prove Proposition [1] below) 
we slightly revise his approach. 

Definition 1 Let b > 2, s > 1, and < u < m be integers and let e = 
(ei, . . . , e s ) £ N s be an s-tuple of positive integers. A (u, m, e, s)-net in base 
b is a point set V of b m points in [0, l) s such that A(J; V) = b m Vo\(J) for 
every elementary interval J in base b of the form 

s 

J = ^M" c^ ^( a ^ + 1 ) rd, ) 

i=l 

with integers di > 0, < < b di , and e^di for 1 < i < s and with 
Vol (J) > b u - m . 

Remark 1 The classical concept of a (u, m, s)-net in base b corresponds 
to the special case e = (1, . . . , 1) in Definition [TJ A simple and well-known 
propagation rule for nets states that a (u, m, s)-net in base b is also a (v, m, s)- 
net in base b for any integer v with u < v < m (see [31 Remark 4.9(2)] and 
[TBI Remark 4.3]). It is trivial that this propagation rule is valid also for 
(u, m, e, s)-nets in base b. 

Remark 2 Tezuka [26] introduced Definition [TJ but he used the condition 
Vol( J) = b u ~ m instead of Vol( J) > b u ~ m . With this original definition, the 
propagation rule for (u, m, e, s)-nets in base b mentioned in Remark [TJ need 
not hold. For instance, choose b = s = 2 and e = (2, 3) £ N 2 . Furthermore, 
fix an integer m > 3. Now take u = m — 3. Then the only possible elementary 
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intervals J C [0, l) 2 in base 2 in Definition [T] with Vol(J) = 2 3 are of the 
form 

J= [0,1) x [a 2 /2 3 ,(a 2 + l)/2 3 ). 

Thus, the net property according to Tezuka's original definition depends only 
on the second coordinates in the point set V. Next take u = m — 2. Then 
the only possible elementary intervals J C [0, l) 2 in base 2 in Definition [T] 
with Vol( J) = 2~ 2 are of the form 

J = [ai/2 2 , (ai + l)/2 2 ) x [0,1). 

Thus, the net property according to Tezuka's original definition depends 
only on the first coordinates in V. It is now clear that with Tezuka's original 
definition, V can be a (u, m, e, 2)-net in base 2 for u = m — 3, but not for 
u = m — 2. 

Definition 2 Let b > 2, s > 1, and w > be integers and let e £ N s . A 
sequence a? , a?!, ... of points in [0, l] s is called a (w, e, s) -sequence in base b if 
for all integers k > and m > u the points [£e n ]& im with fcfe m < n < + l)b m 
form a (m, m, e, s)-net in base 6. 

Just like in Remark [Tj it is clear that the classical concept of a (u, s)- 
sequence in base b corresponds to the case e = (1, . . . , 1) in Definition [2j 

Remark 3 A further advantage of Definitions[T]and[2]is that we can drop the 
additional condition on u and m in Tezuka's original definitions of (u, m, e, s)- 
nets and (u, e, s)-sequences in base b, namely that m — u is of the form 
jiei + ■ — h j s e s with integers j±, . . . , j s > 0, where as above e = (ei, . . . , e s ). 

Proposition 1 Let b > 2, s > 1, and < u < m be integers and let 
e = (ei,...,e s ) GN S . Then any {u, m, e, s)-net in base b is a (t, m, s)-net in 
base b with 

s 

t = min(-u + ^^(ej — 1), m). 

i=i 

Proof. Let V be a (u, m, e, s)-net in base b. If u + J2i=i( e i ~ 1) — m ' 
statement is trivial. Now assume that u + X^=i( e * — 1) < We consider 
an elementary interval J C [0, l) s in base b of the form 

s 

J = l[[a i b- d \{a i + l)b^) 
i=i 



5 



with Vol(J) > 6*" m , i.e., with 

s s 

di < m — t = m — u — ^^(ej — 1). 



i=l i=l 



For each z = 1, . . . , s, let be the unique integer with < Wi < — 1 such 
that g?, + w i is a multiple of e,. Then 



[Oife - *, (a< + = |J [(a^ ! + /i - 1)6-*-*"*, (a^ ! + h)b~ di - Wi ) 



h=l 



is the disjoint union of 6 Wi one-dimensional elementary intervals in base h 
of length b~ di ~ Wi . Consequently, J is the disjoint union of s-dimensional 
elementary intervals Kh, 1 < h < h,Us =: if, in base 6 with 

s 

Vol(K h ) = Y[b~ di ~ Wi = &-£f=i d i-£l=i™i > ft-' 7l + u +EJ=i( e i-i- w i) > 

for I < h < H, where we used ([T]) and Wj < — 1 for 1 < i < s. By 
Definition [TJ we have 

P) = 6 m Vol(i^) for 1 < h < H. 

Therefore 

H H 

A(J;V) = ^A(K h ;V) = 6 m ^Vol(^) = 6 m Vol(J), 

h=l h=l 

and so V is a (£, m, s)-net in base b. □ 

Corollary 1 Let b > 2, s > 1, and u > be integers and lete = (e±, . . . , e s ) G 
N s . Then any (u, e, s) -sequence in base b is a (t, s) -sequence in base b with 

s 

t = u + ^{e i -l). 

i=i 

Remark 4 Tezuka [26] has shown with his original definition of a (it, e, s)- 
sequence that a generalized Niederreiter sequence is a (0, e, s)-sequence in 
base q, where e = (ei,...,e s ) G N s and ej, 1 < z < s, is the degree of 
the ith base polynomial over W q in the construction of generalized Nieder- 
reiter sequences. It is immediately seen from his proof that a generalized 
Niederreiter sequence is also a (0, e, s)-sequence in base q in the sense of 
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our Definition [2J Corollary [T] implies that it is a (t, s)-sequence in base q 
with t = Xli=i( e * — !)• This gives another proof of this well-known result 
from [23]. Consider now the case of the Niederreiter sequences constructed 
in [13] . They form a subfamily of the generalized Niederreiter sequences, and 
so a Niederreiter sequence is also a (0, e, s)-sequence in base q and a (t, s)- 
sequence in base q with t = J2i=i( e i ~ !)• Dick and Niederreiter [2] proved 
that this t- value is the exact value of the quality parameter t for Niederreiter 
sequences, that is, for any Niederreiter sequence no smaller value of t is pos- 
sible. The case of Niederreiter sequences shows therefore that the formula 
for t in Corollary [1] is in general best possible. 

The parameter u of a (u, m, e, s)-net in base q constructed by the digi- 
tal method can be derived from the following result. We use the standard 
convention that an empty set of vectors from a vector space is linearly inde- 
pendent. 

Proposition 2 The matrices C^ 1 ', . . . , E W q nxm generate a (u, m, e, s)- 
net in base q if and only if, for any integers d±,...,d s > with ei\di for 
1 < i < s and d\ + • ■ ■ + d s < m — u, the collection of d\ + ■ ■ • + d s vectors 
obtained by taking the first di row vectors of for 1 < i < s is linearly 
independent over ¥ q . 

Proof. This follows by a straightforward adaptation of standard arguments 
that can be found in [31 Section 4.4.2] and [161 Section 4.3]. □ 

Using Proposition |2j it is easy to determine the parameter u of a (u, e, s)- 
sequence in base q obtained by the digital method. 

Propositions The matrices C^ 1 ' , . . . , C^ s ' E ^ xN o generate a (u,e,s)- 
sequence in base q if and only if, for every integer m > u and all integers 
di, . . . , d s > with Ci\di for 1 < i < s and 1 < d% + • • • + d s < m — u, the 
(di + ■ ■ ■ + d s ) x m matrix over ¥ q formed by the row vectors 

(c {i) c {i) c {i) ) E ¥ m 

with 1 < j < di and 1 < i < s has rank d\ + • • • + d s . 

3 A construction of (u, e, 5)-sequences 

So far, the only verified family of (w, e, s)-sequences with e ^ (1, . . . , 1) is that 
of generalized Niederreiter sequences, according to a result of Tezuka [26] . 
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The construction of generalized Niederreiter sequences in [24], like the earlier 
construction of Niederreiter sequences in [14], is based on rational function 
fields over finite fields. In this section we present a new construction of 
(u, e, s)-sequences using arbitrary global function fields. Recall that a global 
function field F is a finite extension of the rational function field ¥ q (x) for 
some finite field F ? . If ¥ q is algebraically closed in F, then ¥ q is called the 
full constant field of F. 

First we collect some basic facts and notation for global function fields. 
We refer to the monograph [23] for general background on global function 
fields. Let F be a global function field with full constant field ¥ q and genus 
g and let ¥ F denote the set of all places of F. We write deg(P) for the degree 
of a place PgPf and vp for the normalized discrete valuation corresponding 
to P. A divisor D of F is a formal sum 

D= n P p 

PeP F 

with np G Z for all P G ¥p and all but finitely many np = 0. We write also 
np = up(D). The degree deg(D) of a divisor D is defined by 

deg(P) = ^pdeg(P) = M^)deg(P). 

PeP F PeP F 

We say that D x < P 2 if Pp{D x ) < v P {D 2 ) for all P G P F . The principal 
divisor div(/) of / G F* is defined by 

div(/) = Mm 

PeP F 

The Riemann-Roch space 

£(P):={/GP*:div(/) + P>0}U{0} 

is a finite-dimensional vector space over ¥ q . We write ^(P) for the dimension 
of this vector space. Obviously, £(Pi) is a subspace of C(D 2 ) whenever 
Pi < P2- The celebrated Riemann-Roch Theorem [231 Theorem 1.5.15] 
ensures that £(D) > deg(P) + 1 — g, and equality holds whenever deg(P) > 
2g — 1 (see [231 Theorem 1.5.17]). For a rational place P (i.e. deg(P) = 1), 
the Weierstrass Gap Theorem [231 Theorem 1.6.8] says that there are exactly 
g gap numbers 1 = %\ < ■ ■ ■ < i g < 2g — 1, that is, integers ij, 1 < j • < g, 
such that £(ijP) = (-{{ij — 1)P)- Note that for integers n > 2g it is clear by 
the Riemann-Roch Theorem that £(nP) = £((n — 1)P) + 1. 
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Finally, choosing a place P G Fp with degree e and a local parameter z 
at P, that is, an element z G F such that z/p(z) = 1, then any / G F has a 
/oca/ expansion at the place P of the form 

oo 

/ = E 

fc=feo 

where /c G Z with Vp(f) > k and G F ?e for all k > k . A detailed 
description of how to obtain a local expansion can be found in [201 pp. 5-6]. 

For our construction we choose a dimension s G N, a finite field W q , a 
global function field F with full constant field ¥ q and genus g, and s + 1 
distinct places P^, Pi, . . . , P s of F satisfying deg(Poo) = 1. The degrees of 
the places Pi, . . . , P s are arbitrary and we put = deg(Pj) for 1 < i < s. 
First we construct a sequence (y r ) r >o of elements of P using the information 
on the dimensions £(nP) for n G No collected above. This allows us to choose 
for each integer r > an element 

y r G £(n r Poo) \ C((n r - l)Poo)- 

Here no < ni < n2 < • • • are the elements in increasing order of No \ 
{«!,..., i g }, where zi, . . . , i g are the g gap numbers of P^. Note that n r = 
r + g for r > g. From the construction of the y r we see that vp(y r ) > for all 
places P G P_f \ {Poo} an d all r > 0. Thus for 1 < i < s, the local expansion 
of y r at Pj (with local parameter at Pj) has the form 

oo 

^ = E^" with all /3«GF^. 

fc=0 

Using an ordered basis of /F 9 , we can identify (3^ r with a column vector 
6? G F e \ 

Now we construct generating matrices 1 < i < s, over ¥ q colum- 
nwise by defining column 0, column 1, etc. For column r (r > 0) of 
we concatenate the column vectors b^ r , k = 0, 1, . . .. This completes the 
construction. 

Theorem 1 Let F be a global function field with full constant field ¥ q and 
let Poo, Pi, . . . , P s be s + 1 distinct places of F with deg(Poo) = 1. Then 
the matrices . . . , G F^ xNo constructed above generate a [u,e,s)- 
sequence in base q with u = g and e = (e±, . . . ,e s ), where g is the genus of 
F and = deg(p) for 1 < i < s. 
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Proof. We write C® = (cfl) £ F^ xN ° for 1 < i < s. According 

^ ' 3>l,r>Q 

to Proposition [31 it suffices to show that for every integer m > g and all 
di, . . . , d s 6 M such that ei\di, . . . , e s \d s and 1 < di + • — h d s < m — g, the 
(di + • • • + d s ) x m matrix M over ¥ q formed by the row vectors 

with 1 < j < di and 1 < 2 < s has rank c?i + ■ ■ • + d s . 

In order to prove this assertion, we show that the kernel of M has dimen- 
sion <m—(di + --- + d s ). We take an arbitrary element (v , • • • , ^m-i) £ F™ 
of the kernel of M. Then 

m—1 

UrC W =0 for 1 < j < d u 1 < % < s. (2) 



■r=0 



In the following, we make use of the element / := YlT=o Vr Vr °f F. From the 
construction of the y r and the definition of / we derive that 

f E £ (n m -iPx>) , 

which is a Riemann-Roch space over ¥ q with dimension m since the divisor 
n m _xPoo = (m + g — i)Poo has degree m + g — 1 > 2g — 1. Note that 
y , . . . ,y m -i form a basis of this Riemann-Roch space, and so each element 
of £ (n m -iPoo) has a unique representation as an F^-linear combination of 
2/o, • • • ,2/m-l- 

From (j2J) and the construction of the matrices . . . , C^ 8 > we deduce 
that 

v Pi {f)>— forl<i<s. 



Therefore 



.s 



which is a subspace of £^n m _iPooj. The degree of the divisor 

-D := ™m-lPx> - ^ — -P 

i=l e * 

satisfies 

deg(P) = m + 5- - 1 - {d x H h d a ) > 2g - 1. 
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By the Riemann-Roch Theorem we know that the dimension of C(D) is 
m—(di + - ■ -+d s ). In view of the one-to-one correspondence between elements 
of F™ and elements of C (n m _iPoo) noted above, this yields the desired upper 
bound on the dimension of the kernel of M. □ 

Corollary 2 With the conditions and notation in Theorem [21 the matri- 
ces C^ l \ . . . , C^ s ' constructed above generate a (t, s) -sequence in base q with 
quality parameter t = g + J2i=i( e i ~~ !)■ 

Proof. This follows immediately from Theorem [1] by applying Corollary [TJ □ 

Example 1 If F — ¥ q (x) is the rational function field over ¥ q (here g = 
0) and Pqo is the infinite place of F, then the sequence (y r ) r >o consists of 
elements of the polynomial ring ¥ q [x] and we have deg(y r ) = r for all r > 0. 
Furthermore, if Pi,...,P s are places of F corresponding to distinct monic 
irreducible polynomials over ¥ q (see [23j Section 1.2]), which serve as local 
parameters in the local expansions, then our construction is consistent with 
the one introduced in [5]. 

Remark 5 So far, the known constructions of generating matrices for digital 
(t, s)-sequences using arbitrary global function fields proceeded rowwise. The 
rowwise approach constructs for each component i G {l,...,s}a sequence of 
elements in the global function field via Riemann-Roch spaces [TBI EZ] or via 
differentials [T2]. Then the row entries of the generating matrices are deter- 
mined by using the coefficients of the local expansions of these elements at 
a fixed rational place. Our columnwise method constructs a sequence of ele- 
ments in the global function field using Riemann-Roch spaces that are built 
using a fixed rational place. Then the columns of the ith generating matrix 
are determined by using the coefficients of the local expansions at the 
place Pi corresponding to the ith component. The first columnwise construc- 
tion of generating matrices was introduced by Niederreiter |15j . which was 
later generalized by Hofer 18] . Those constructions exploited properties of 
certain elements of the polynomial ring F g [a;]. 

Remark 6 The basic ingredients of the construction of Niederreiter-Xing 
sequences in [27] are, as in our construction, a global function field F with full 
constant field ¥ q and s+1 distinct places P^, P±, . . . , P s of F with deg(Poo) = 
1 and the degrees = deg(Pj) for 1 < i < s being arbitrary. We claim that 
a Niederreiter-Xing sequence with these data is a (g, e, s)-sequence in base 
q, where g is the genus of F and e = [e\, . . . ,e s ). This is shown by an 
appropriate modification of the proof of [2JJ Theorem 1]. We adhere to the 
notation in that proof. In view of Proposition [3J it suffices to show that 
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for every integer m > g and all d±, . . . , d s € No such that ei\di, . . . , e s \d s 
and 1 < di + ■ ■ ■ + d s < m — g, the vectors 7r m (c^) G F™, 1 < j < d i} 
1 <i < s, are linearly independent over ¥ q . In the course of the proof of [2"7l 
Theorem 1], an auxiliary element 

M» + E([^H p ')= £ ( D+ i:t*) 

is considered. Furthermore, it is shown that vp^ik) > m + g + 1. These two 
facts, together with vp^D) = 0, imply that 



k e C(D + ~ p i - ( m + 9 + !) p o. 

i=l * 



Since deg(-D) = 2g, we have 



det 



g ( D + J2-Pi-(m + g + 1)P 



1 < 0, 



8=1 



and so k = by [23, Corollary 1.4.12(b)]. An appeal to [271 Lemma 2] 
completes the proof. By Corollary [TJ the Niederreiter-Xing sequence under 
consideration is a (i, s)-sequence in base q with i = g + J2i=i( e i ~~ This 
recovers the main result of [27] in a different way. This value of t is the 
same as in Corollary [2j However, it should be pointed out that our new 
construction is considerably simpler than the one in [27J. For instance, the 
auxiliary positive divisor D of F with deg(-D) = 2g and vp x (D) = in 
is not needed in our construction. 



4 Finite-row generating matrices for (u, e, s)- 
sequences 

Columnwise constructions as in this paper have certain advantages, for in- 
stance, one may obtain finite-row generating matrices more easily with this 
approach. For more details and the motivation for finite-row generating ma- 
trices, we refer the interested reader to [7J, [8] and the references therein. To 
drive home an important point, we cannot refrain from emphasizing one great 
asset of finite-row generating matrices, namely that the sparsity of finite-row 
generating matrices leads to a considerable speedup in the computation of 
the matrix-vector products that form the main operation in the actual im- 
plementation of the digital method (see Section [p. 
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In this section we explore a more refined version of the columnwise con- 
struction in Section [3] using global function fields. We choose the y r in the 
sequence (y r ) r >o by using subspaces of £(n r Poo) and in this way we are able 
to construct finite-row generating matrices with, in a certain sense, optimal 
row lengths. 

As in Section [3J let s G N be a given dimension, let ¥ q be a finite field, 
let P be a global function field with full constant field ¥ q and genus g, and 
let Pqo, Pi, . . . , P s be s + 1 distinct places of F satisfying deg(Poo) = 1 and 
deg(p) = e, for 1 < i < s. We construct (y r ) r >o as before, but in the case 
where r > g we choose y r such that 



Here v := lcm(ei, . . . , e s ), whereas the integers U > are determined by 
division with remainder in the nonnegative integers as follows: 

r - g = svl s + w s , w s G {0, . . . , sv - 1}, 

w s = (s - l)vl s -i + W 8 -t, W s -i G {0, . . . , (s — l)v — 1}, 

W 2 = Vlx + Wi, Wi G {0, . . . ,v — 1}. 

Note that such a y r exists since D' r < D r and deg(P r ) > r + g — (r — g) = 2g 
and deg(P',) = deg(P r ) — 1. The rest of the construction is carried out as in 
Section [3] and leads to the new generating matrices . . . , G F^ xN °. 

From the proof of Theorem [1] and from Corollary [2j we see that the 
generating matrices . . . , constructed in this way within the digital 
method yield a (g, e, s)-sequence in base q with e = (ei, . . . , e s ) and also a 
(t, s)-sequence in base q with t = g + J2i=i( e i ~ !)■ O ur more refined choice 
of the sequence (y r ) r >o of elements of P leads to bounds on the number of 
nonzero entries in each row of the generating matrices. To describe these 
bounds, we define the length of a row (co,ci, . . .) G F^ with only finitely 
many nonzero entries by max {7 G N : Cj_i 7^ 0}, with the proviso that the 
length is if c r = for all r > 0. 

Theorem 2 The generating matrices . . . , of the (g, e, s) -sequence 
in base q constructed in this section satisfy the property that for all 1 < i < s 




and 
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and (ieN, the length of the dth row of is at most 



d-1 



g + sv 



+ iv. 



v 



Here g is the genus of F , e — (ei, . . . , e s ) with e { = deg(Pj) for 1 < i < s, 
and v = lcm(ei , . . . , e s ) . 

Proof We fix i G {1, . . . , s} and d G N and let a G N and w G {1, . . . , v} 
be such that d = av + w. We write C w = (cf r ) G F^ xN °. In order 

to ensure that for every r > g + st> |_(d — 1)A"J + iv = g + sva + if, the 
entry c^ r in the dth row of is 0, it suffices to show by the construction of 
that up-(y r ) > (a + l)v/ei. (Note that (a + l)u > d.) Preliminarily, we 
consider the case where r — g > sva and we see from the definition of l s that 
Z s > a. Now in the case where r — g > sva + iv, we have either i s > a + 1, 
or otherwise l s = a and £j > 1 for at least one j G {«, i + 1, . . . , s — 1}. Thus, 
we always have k + + • ■ ■ + l s > a + 1. From w r G C(D r ) we deduce that 



Remark 7 In the case where the genus g is and t> = 1, the matrices 
. . . , constructed in this section have asymptotically shortest possi- 
ble row lengths in the sense of Hofer and Larcher [91 Section 4] (compare also 
with p. 590], Remark 3], and [II]). 

Remark 8 Given a prime power q and a dimension s G N, let F be a global 
function field with full constant field F g such that F contains at least s + 1 
rational places. Then in the construction in this section we can, besides P^, 
also take P\, . . . , P s to be rational places of F, that is, = deg(Pj) = 1 for 
1 < i < s. In this way we obtain a (t, s)-sequence in base q with finite-row 
generating matrices and with t = g, the genus of F. Now we optimize the 
construction in the following sense. As in [TSJ Eq. (10)], we define V q (s) to 
be the least g G No such that there exists a global function field with full 
constant field W q , genus g, and at least s + 1 rational places. Then there 
exists a (V^(s), s)-sequence in base q with finite-row generating matrices. It 
was shown in [T51 Theorem 4] that V g (s) = O(s) as s — > oo with an absolute 
implied constant. Hence for fixed q, we obtain for any dimension s a (t, s)- 
sequence in base q with finite-row generating matrices such that the value 
of the quality parameter t grows linearly in s as s — > oo. This growth rate 



VPiiVr) > {k + k+i H h h)~ > 



(a + l)v 



and this proves our claim. 



□ 
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of t is best possible since it is well known that for any fixed base 6, the 
least t- value t b (s) of any (£, s)-sequence in base b grows at least linearly as 
a function of s as s — > oo (see [Hi Theorem 8] and [21,, Theorem 1]). For 
the only construction of (t, s)-sequences in base q with finite-row generating 
matrices available for any q and s that was known previously, namely the 
construction in [S], the least possible value of t grows at the rate slogs as 
s — > oo for fixed q. 



5 Discrepancy bounds 

All known upper bounds on the star discrepancy D^(S) of a (t, s)-sequence 
S in base b are of the form 

D* N {S) < cN-^logN) 3 + O {N- 1 (log N) 3 - 1 ) for all N>2, (3) 

where the constant c > and the implied constant in the Landau symbol 
depend only on b, s, and t. The classical formulas for c were established by 
Niederreiter [131 Section 4] and they are summarized in [IBJ Theorem 4.17]. 
Later, ameliorated values of c were obtained by Kritzer [10J and Faure and 
Lemieux [6]. The currently best values of c are those of Faure and Kritzer [5], 
namely c = cfk given by 



cfk 




if b is even, 
if b is odd. 



It follows that for every base b > 2 we have 



b l 1/6—1 . 

CFK > -7 ■ 'A 7TT-T • (4) 



s! 2 V.21og6 

By applying the signed-splitting technique of Atanassov [1], Tezuka 
established an upper bound on the star discrepancy D* N {$) of a (u, e, s)- 
sequence S in base b which is of the form (j2J) with the constant c = CTez given 
by 

c _ hU TT L &6 V 2 J < & " b ei+ - +e ° ^ 
s ' e « 1°S ^ _ s ' (2 log6) s ei • • • e s 

Note that our Definition [2] of a (u, e, s)-sequence in base b is slightly stronger 
than Tezuka's original definition in [2B], and so Tezuka's discrepancy bound 
is obviously valid for our concept of a (u, e, s)-sequence in base b as well. 

The new sequences constructed in Sections [3] and H] are (t, s)-sequences 
in base q as well as (u, e, s)-sequences in base q with suitable t, u, and e. 
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Therefore we can apply two versions of the discrepancy bound ((3]) , namely the 
one with c = cfk and the one with c = c^ ez . It is of interest to compare these 
two bounds. Recall that for these sequences we have u = g, the genus of the 
global function field F, as well as e = (e±, . . . , e s ) and t = g + J2i=i( e i ~ 1)> 
where e« = deg(Pj) for 1 < % < s. Then from (J4]) and ([5]) we obtain 

c FK 1 (q - 1 



c Tez \ i / . 



^mn- 

Thus, if 



then the discrepancy bound with the constant CTez is better than the one 
with the constant cfk- The condition ([6]) will be satisfied in many cases. 

Consider a typical situation where we fix a global function field F with 
full constant field ¥ q such that F has at least one rational place P^. We 
arrange all distinct places in Fp \ {-Poo} into a list Pi, P2, . . . in an arbitrary 
manner. Then we have the following result. 

Proposition 4 Let P\,P2, ... be a list of all distinct places in Fp \ {Poo} 
as above, where F is a given global function field with full constant field ¥ q . 
Then for any < e < 1/q we have 



ndeg(P)>(log gS )« 1 ^- ^ 



for all sufficiently large s, where log q denotes the logarithm to the base q. 

Proof. For any k G N, let P>k be the number of places of F of degree k and 
let Tk = Ylh=i Bh be the number of places of F of degree < k. Let fo = 
deg(Poo) = l,/x,/2, ... be the nondecreasing sequence of positive integers 
that is obtained by listing each k G N with multiplicity Note that 
fo, fii fi-, ■ ■ ■ is also the list of the degrees of all places of F in nondecreasing 
order. 

Throughout the proof we assume that s is sufficiently large. Let j = j(s) 
be the largest h G N with < s + 1. Then 

s s s j 

log deg(P) > ]T log fi = fi>Y, log fi = E Bk lo § k - ( ? ) 

i=\ i=l i=0 i=0 k=l 
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For any < 5 < 1 we have 



B k log k > B * lo § k ^ 5{T i - T m-J 
k=i k =[j s ] 

The Prime Number Theorem for global function fields pj] yields 

q qi ( qi 



Tj = • ho — ) as j — > oo. (9) 

3 q-1 j \ j ) 

Noting that with s also j = j(s) is sufficiently large, we deduce from flU), 
and (jHJ) that for any < S\ < q/(q — 1) we have 



a j / 

]T log deg(P) > B k log k > I — 

i=l k=l ^ 



The definition of j implies that T J+1 > s + 1. Using again (jH]), we obtain for 
any e 2 > that 

, <T , +I <f_L +ei ^<f_3L +ei ,W 



g - l y j + l \g - l 

and so for any < e 3 < (q — l)/q 2 we have 
Similarly, 



s < T j+1 < — — < q\ 
J + 1 

and so j > \og q s. By combining these inequalities, we get 

Y lo S deg(Pi) > ^^TT _ £l ) \ ~ £s ) S log log 9 S ' 

Now for a given < e < 1/q we can choose S\ and £3 suitably small such 
that 

' S—-e\ ( q ~ l - s 3 ) > - - e 
o-l V V 9 2 7~ 9 



Consequently, 



^logdeg(P) > Q -ej sloglog 9 s 



for all sufficiently large s, and this implies the desired result. □ 
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If we now take for the given global function field F and a given s G N 
the first s places Pi, . . . ,P S from the list Pi, P2, . . . in Proposition HJ then it 
follows from this result that the condition ([6]) is satisfied for all sufficiently 
large s. 

Remark 9 In view of Remark [6j the discussion in this section applies in ex- 
actly the same way to upper bounds on the star discrepancy of the Niederreiter- 
Xing sequences constructed in [27] . 
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